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Abstract 

For a two-dimensional system of agents modeled by molecular dynamics, we si- 
mulate epidemics spreading, which was recently studied on complex networks. Our 
resulting network model is time-evolving. We study the transitions to spreading as 
function of density, temperature and infection time. In addition, we analyze the 
epidemic threshold associated to a power-law distribution of infection times. 
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1 Introduction 



The statistical spreading of infections, information or damage, involves non- 
equilibrium phenomena. Fluctuations and spatial correlations play an impor- 
tant role and are often not solvable exactly. Usually these processes are studied 
on a lattice that can be regular [1,2], hierarchical or small- world [3]. But in 
most cases the population in question is mobile. Therefore, in this work we 
study a system of particles moving according to a simple Newtonian dynam- 
ics. We simulate on it a known contact process [4], described in terms of a 
'SIS' model of infection, or infection without immunization, that is the state 
of the particles are healthy or infected, and are susceptible to re-infection after 
healing, thus the name of the model (SIS: susceptible-infected-susceptible). 
We characterize the transition to spreading of the epidemic dynamics, and 
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obtain a continuous range of critical exponents changing the density of the 
system. The observed behavior results to be the 'SIS' analogous of a model 
of 'stirred percolation' [6,5,7] which was used to describe epidemic dynamics 
with immunization ('SIR'). 



2 Model 



The simulations are carried out on a square shaped cell of linear size L with 
periodic boundary conditions. L = ^Jp/N, is given by the number of particles 
(N) and the density (p). The particles are represented by 'soft-disks' of radius 
ro moving continuously on the plane. The interaction between two particles 
at positions r*j and Vj corresponds to a Lennard- Jones potential truncated at 
its minimum, 

+ U Q , \ri- rj \ < 2r ,(l) 
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reduced units are used in which U , r , k B (the Boltzmann constant) and m 
(the particle mass) are all unity. 

Along this paper, the particles are considered to be 'agents', to model the 
'SIS' process described above. In a simple version one can assume that each 
time two agents interact or 'collide' (that is, if \ri — Vj\ < 2r ), the 

infection propagates from an infected agent to a susceptible one. In most of 
our simulations we used a simple initial condition: (i) at time t — 0, N agents 
are distributed regularly in the cell, (ii) have the same absolute velocity v with 
randomly distributed directions, (iii) a central agent is infected and the rest 
are susceptible. At each time step the positions {r ; }, velocities {v ; } and the 
infection state {a,j} of the system are updated. We use the molecular dynamics 
(MD) scheme of cell subdivision with the leapfrog integration method [8]. 
Once an agent is infected, it heals and becomes susceptible again after a fixed 
number of time steps, that is the 'time of the infection' (At in f), and is a free 
parameter of the model. 



3 Scaling 



The dynamics of an epidemics is described in terms of the infection rate (A). 
That is defined as the number of agents one agent infects before healing. 
Therefore, in this model 

A = At inf /t co u, (2) 
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Fig. 1. Starting with one infected agent, the figures shows how the infection spreads 
with time and after a certain period of time ('transient'), the fraction of infected 
agents (Fjm) fluctuates around a mean value. The fraction of infected agents in this 
'quasi-equilibrium' state, increases with the infection rate (A). 

where t co u is the mean time between two collisions, and depends on the mean 
velocity of agents {{v)). Neglecting the interaction potential with respect to 
the kinetic energy, one has: 



(3) 




Thus, the mean number of collisions {{n coU )) of one agent during a period of 
time t, is given by the area within which it interacts (2r (v )t), multiplied by 
the density, 

(n coU ) = p 2r (v)t. (4) 



The pre-f actor of time in eq. (4) is the collision frequency, and its inverse, 
gives the mean time between collisions, 



Tcoll ~ P 2rJjrk B (5) 



In Fig. 1 we see the fraction of infected agents (FjM{t) vs. time) for different 
values of temperature (T), density (p) and time of infection (At in f). In each 
case we start with one infected agent and after a transient, the system fluc- 
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Fig. 2. The plot shows the collapse of data of Fig. 1. The vertical axis was divided 
by the mean field approximation of FjmW and the horizontal axis is divided by 
the collision time (t co u). 



tuates around a value of Fjm, which depends on A. This mean value can be 
calculated using a mean field approach, 

= -F IM (t) + XF IM [1 - F IM ] (6) 



The first term of the r.h.s is the fraction of agents that heals and the second 
term, is the fraction of agents that becomes infected. After the transient, 
dF IM (t)/dt ~ 0. Thus 



,0 if A < A c 

%W= ~ ( 7 ) 

1 - 1/A if A > A c 



where A c = 1, is known as the critical rate of infection. In Fig. 2 we see the 
collapse of the realizations of Fig. 1, with the expression obtained in eq. (7). 
The dynamics of the system is characterized by the infection rate A, which 
contains all the free parameters of the system. 
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Fig. 3. Fraction of infected nodes Fjm as function of A for different conditions 
of T and p. The data agree with the MF approximation (eq. (7)) only for low 
densities (< 0.1). The inset shows the log — log plot of Fim vs. A — A c with T = 1 
and increasing densities: p = 0.05, A c = 1.05 (diamonds), p = 0.1, A c = 1.06 
(circles), p = 0.23, A c = 0.862 (filled triangles), p = 0.30, A c = 0.73 (crosses), and 
p = 0.46, A c = 1.05 (rotated triangles). The lines are regressions of the data and 
have respectively slopes: = 1.0, = 0.92, = 0.738, = 0.66, and = 0.599. 

4 Transition to spreading 



We analyze the transition to spreading of the epidemic. In systems of 1254 
agents Fjm after the transient is averaged, over 500 initial conditions at a 
given A. In Fig. 3 we see that the results agree with the mean field approxi- 
mation for densities lower than 0.1. Increasing the density, changes the shape 
of the transition curve. Near the critical point the Fjm follows a power law, 
Fim ~ (A — Ac)' 3 . The inset of Fig. 3 in a double-logarithmic plot, demonstrates 
this power law behavior, with straight lines of slope (3. As can be seen, the 
value of P depends on the density of the system. Smaller values of (3 indicate 
more significant changes of -FVm(A) near the transition. Even for low densities 
(i.e. p G [0.1,0.2]), where one expects MF to be valid, the critical exponents 
change with the density. 

In order to confirm this observation, we study the transition in detail. Fig. 4 
shows the log — log of Fjm vs. (A — A c ), averaging over 1000 realizations in 
systems with p = 0.2, T = 1, and A c = 0.935. We see that the data fit the 
expression ~ (A — \ C Y , with exponent = 0.773. 
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Fig. 4. Log — log plot of the fraction of infected nodes Fjm as function of A — A c 
for (p = 0.2, T = 1), with A c = 0.935. The dashed line is a fit to the form 
Fim ~ (A - Ac)^, with exponent (3 = 0.773. 
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Fig. 5. Evolution of the survival probability of infection P(t), the mean number of 
infected agents n(t) and the mean square distance of spreading R 2 {t) in time. Each 
graph contains three curves near criticality. With p = 0.2, from bottom to top: 
A = 0.93, 0.94 and 0.95. The circles come from regressions to calculate the critical 
exponents 

Other critical exponents, are obtained, if one averages the survival probability 
of infection P(t), the number of infected agents n(t) and the square distance 
of spreading R 2 (t) and plots them against time. At the critical point, they are 
expected to display asymptotic power laws [9], 



P(t) ~ t~ 5 , n(t) ~ t\ R 2 {t) ~ t z . 



(8) 
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MF ,0 = 0.1 p = 0.2 Contact Process(2D) [2] 
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1.0(8) 
0.9(2) 
0.5(9) 
0.1(5) 
1.3(0) 



0.9(4) 
0.7(7) 
0.5(3) 
0.2(5) 
1.2(7) 



1.6488(1) 
0.583(4) 
0.4505(10) 
0.2295(10) 
1.1325(10) 



Table 1 

Critical rate of spreading and exponents for SIS model on moving agents with 
p = 0.1 and 0.2, contact process in two dimensions, and estimates obtained by 
mean field MF. 



1.00 



0.90 



0.80 



0.70 



0.60 



Mean Field 



Contact Process 2D 



0.50 ' — 
0.00 



0.10 



0.20 



0.30 



0.40 



0.50 



Fig. 6. Numerical estimates for the critical exponent 0, for systems with different 
density. 

The relations in eq. (8) apply at long times, and require that the infection 
does not reach the boundaries of the system. Results for the three quantities 
P(t), n(t) and R 2 (t), averaged over ~ 10 3 realizations, with systems of ~ 10 4 
agents, and fixed temperature (T — 1), are shown in Fig. 5 for p = 0.2. 

The values of the critical exponents are reported in table 1. The known hyper- 
scaling relation of the dynamics 45 + 2i] = dz, where d is the dimension of 
the system, is recovered within the range of numerical errors. Increasing the 
density of the system gives a continuous change in the critical exponents of 
the epidemic dynamics, they go from MF values to the exponents for contact 
process in a two dimensional lattice [2]. Numerical estimates of the critical 
exponent j3 vs. p are shown in Fig. 6. 



The dependence of the critical exponents of the dynamics on the density of 
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agents (p) , is analogous to the dependence on the 'flammable' fraction of space 
((f)), observed in 'forest fire' dynamics, or epidemics dynamics with immuniza- 
tion [5,7], described through 'stirred percolation' models. 'Stirred percolation' 
consists in random walkers on a changing lattice and was proposed by de 
Gennes [6] to describe the behavior of conductivity in binary mixtures [10]. 



5 Power-law distribution of infection times 

Considering the same value of infection times (At in f) for each agent corres- 
ponds to situations with homogeneous connectivity. In order to extent the 
model to some real world situations where the number of contacts varies 
greatly from one agent to another, we assign At in f to each agent following 
a power-law distribution, 

P(At inf ) = ( 7 - 1)AQ At mf > 1. (9) 

As a result, for 2 < 7 < 3, the epidemic threshold tends to zero, like has been 
observed for SF networks [3] (see Fig. 7). 

However, the shape of the spreading curve has a point of inflection (p c ) above 
which, the infection is much larger (Fim ~ (p — Pc) 13 , Pc = 0.06(5), (3 = 2.65, 
for 7 = 2.4). 

We see that the more 'connected' agents are responsible for the absence of an 
epidemic threshold (tail of the curve). The infection would spreads only if the 
mean rate of spreading is larger than one, that is 

(At inf ) 7-1 /ttT/cb f . 

A c = — = p 2r \ > 1. 10 

Tcoii 7-2 f m 

A c = 1 gives the inflection point p c , which for 7 = 2.4 is p c = 0.0718 (according 
to eq. 10, for r = 2 1 / 6 and T — 1), and agrees with the numerical value 
p c = 0.06(5) reported in Fig. 7 averaging over about 10 3 realizations with 
about 10 3 agents each. 



6 Conclusions 

Novel effects are observed studying the SIS model of infection on a system 
of moving agents. A continuous range of critical exponents are observed as 
function of the density of agents, recovering mean field predictions for lower 
densities, two-dimensional exponents of contact process, increasing the density. 
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Fig. 7. Fraction of infected nodes vs. density, with power-law distribution of infection 
time (eq. 9) for 7 = 2.4. Fjm ~ {p — Pc) 13 is shown in the upper-left-corner inset 
with p c = 0.065 and (3 = 2.6(5). The bottom-right-corner inset is the main plot in 
log — log scale to see better the tail of the spreading curve, for p < p c . 

Introducing a power-law distribution of infection times, the epidemic threshold 
becomes zero due to the more 'infecting agents'; but still there is a criticalr&te 
of spreading, which depends on the exponent of the distribution and the mean 
time of interaction among the agents. 
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